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Large scale structure surveys will likely become the next leading cosmological probe. In our universe, 
matter perturbations are large on short distances and small at long scales, i.e. strongly coupled in 
the UV and weakly coupled in the IR. To make precise analytical predictions on large scales, we 
develop an effective field theory formulated in terms of an IR effective fluid characterized by several 
parameters, such as speed of sound and viscosity. These parameters, determined by the UV physics 
described by the Boltzmann equation, are measured from A^-body simulations. Wc find that the 
speed of sound of the effective fluid is ~ 10~^c^ and that the viscosity contributions arc of the 
same order. The fluid describes all the relevant physics at long scales k and permits a manifestly 
convergent perturbative expansion in the size of the matter perturbations d{k) for all the observables. 
As an example, we calculate the correction to the power spectrum at order S{k)^. The predictions of 
the effective field theory are found to be in much better agreement with observation than standard 
cosmological perturbation theory, already reaching percent precision at this order up to a relatively 
short scale k ~ 0.24/i Mpc~^. 



1 Introduction 

Large Scale Structure Surveys have the potential of becoming the leading cosmological ob- 
servable in the next decade. They contain a tremendous amount of cosmological information. 
If we were able to extract information from all the modes that go from the horizon scale 
~ 10^ Mpc to the non-linear scale ~ 10 Mpc, we would obtain about 
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independent modes. The Planck satellite in comparison has about (2 x 10^)^ ~ 10^ modes. Of 
course, accessing all this information is much harder than for the CMB, due to the short scale 
non-linearities. There are several aspects to this problem. The first problem is related to our 
currently limited understanding of the evolution of dark matter on large scales. Non-linear 
corrections are very important even on scales larger than 10 Mpc, because modes of different 
wavelengths couple to each other. Understanding these corrections is a problem that affects 
all large scale structure observables. There are then two additional issues that affect most, but 
at least not all, observables. One is the fact that most dark matter is clumped in very non- 
linear structures (dark matter halos); and the other is the fact that what we often observe 
are galaxies, and not just dark matter halos, and not even dark matter long wavelength 
perturbations. The solution to these two last problems requires the correct understanding 
of the so-called halo- and galaxy- biases. These two problems, while important and deeply 
interesting in their own right, are very astrophysical in nature, and we do not address them 
here. 

Instead here we try to address in a rigorous way the first problem, that is the prediction 
of the dark matter distribution on scales larger than the non-linear scale. The fact that the 
universe is characterized by two well separated scales, the Hubble scale, over which pertur- 
bations are linear, and the non-linear scale, which indeed characterizes the scale over which 
gravitational collapse overtakes the expansion of the universe, makes the problem amenable 
to an Effective Field Theory (EFT) treatment. An effective theory is a description of a system 
that captures all the relevant degrees of freedom and describes all the relevant physics at a 
macroscopic scale of interest. The short distance (so called ultraviolet or 'UV') physics is 
integrated out and affects the effective field theory only through various couplings in a per- 
turbative expansion in the ratio of microphysical UV scale/s to the macroscopic scale being 
probed. This technique has been systematically used in particle physics and condensed matter 
physics for many years, but has not been fully used in astrophysics and cosmology. An impor- 
tant early (and recent) application of these techniques in cosmology is the so-called Effective 
Field Theory of Inflation [1]. In a similar vein, understanding the large scale properties of 
the universe is very important, and is ready for a careful analysis. 

Indeed the situation in the universe is very similar to what happens in the chiral La- 
grangian that describes pion interactions in Particle Physics. At very low energies, pions are 
weakly interacting. These interactions and the size of the fluctuations grow with energy until 
we hit the Quantum ChromoDynamics (QCD) scale, ~ 47rF^, at which the pions become 
strongly coupled. The Chiral Lagrangian [2] offers the correct effective theory allowing ar- 
bitrarily precise predictions, up to non-perturbative effects, at energies E <C iirFT^. In our 
universe, matter fluctuations are small at large distances and becomes larger and larger as 
we move up to the non-linear scale. Since the size of the non-linear terms, which are nothing 
but interactions, grows with the size of the fluctuations, we see that at long distances the 
universe should be described by some weakly coupled degree of freedom, that becomes more 
and more interacting as we move closer to the non-linear scale, at which point the fluctuations 
become strongly coupled. The coupling constant should indeed be represented by the ratio 
of the considered wavenumber k over the wavenumber at the non-linear scale k^L: k/k^L. 
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Notice that indeed the size of the density perturbations 5p/ p on a scale k scales as {k/k^LY. 
This scahng suggests the existence of an effective field theory that should allow us to describe 
with arbitrary precision the universe on scales k ^ fcjvL, very much as the Chiral Lagrangian 
represents the right effective field theory to describe pion dynamics to arbitrary precision. 

Such an effective theory would have particularly relevant observational implications. Al- 
ready now, large scale structure surveys such as BOSS or DES are measuring the galaxy- 
galaxy correlation function, so called Baryon Acoustic Oscillations (BAO), at scales of order 
100 Mpc. Next generation experiments such as LSST will measure this quantity at about 
percent precision. These observations contain huge amount of information on Dark Energy 
and on Inflation, through for example the non-Gaussianity of the primordial perturbations. 
The BAO scale is about one order of magnitude longer than the non-linear scale, where 
5p/ p ~ 10^^, and therefore physics at this scale must be describable by rigorous pcrturbativc 
methods. The alternative is to rely on either time consuming numerical simulations, or on 
analytical approaches that however are limited by some irreducible mistake that is hard to 
quantify precisely. In an ideal situation, numerical A'"-body simulations should be quickly 
done only at small scales, to describe phenomena affected by gravitational collapse, rather 
than running large simulations to describe weakly coupled physics. This has been indeed 
recently elucidated in the context of the bias, where it was shown that in order to derive the 
bias on large scales one needs to run very small simulations in a curved universe [3]. This 
line of reasoning is indeed very similar to what happens in QCD, where we perform lattice 
simulation to measure quantities relevant at energies above around one GeV, while we use 
the chiral Lagrangian for predictions at smaller energies. 

The effective field theory (EFT) of the long distance universe was initially developed by 
some of us in [4]. It was noticed that by concentrating on length scales longer than the 
non-linear scale, the universe is described by a fiuid with small perturbations. The equations 
of motion of this fiuid arc organized in a derivative expansion in the ratio of the considered 
wavenumber over the wavenumber associated to the non-linear scale kjsiL ~ 1/10 Mpc~^. At 
leading order in derivatives, the fluid has the stress tensor of an ordinary imperfect fluid, 
characterized by a speed of sound for the fluctuations, a bulk and a shear viscosity, plus 
a stochastic pressure component. This makes our approach different with respect to the 
'standard' approaches both at a quantitative and a qualitative level. 

The purpose of this paper is to further develop this effective theory and be able to make 
observational predictions. The parameters that characterize the fluid, the speed of sounds, 
bulk viscosity, etc., are determined by the microphysics at the non-linear scale, that we call 
UV, and cannot be derived from within the effective theory. They have to be either flt to 
observations, or measured in small A^-body simulations. At this point, the EFT becomes 
predictive. Again, this is very similar to what happens in QCD, where one can measure the 
pion coupling constant F^^ in lattice simulations, after which the Chiral Lagrangian becomes 
predictive. 

Our basic method and key results are summarized as the following: 

• By smoothing the coUisionlcss Boltzmann equation for non-relativist ic matter in an 
expanding FRW background on a length scale A~^, we establish the continuity and 
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Euler equations for an effective fluid. Tfie Euler equation includes an effective stress- 
tensor [t'-'Ja that is sourced by the short- modes Sg. 

• By taking correlation functions of the stress tensor in the presence of long wavelength 
fluctuations, we define an effective stress-tensor that is only a function of the long 
wavelength fluctuations. It takes the form 



where the various parameters c^, cf^ etc. are defined by proper correlation functions of 
short wavelength and long wavelength fluctuations. 

• By directly evaluating the stress tensor from the microphysical theory, i.e., from A^- 
body simulations, and computing the appropriate correlation functions, we calculate 
the value of the fluid parameters. For a ACDM universe with standard cosmological 
parameters at redshift z — and smoothing scale A = l/3h Mpc"^, we find 



where c^^^^}^ is the combination of c^, c^^ and c^^ that is relevant for the leading non- 
linear correction to the power spectrum (one-loop in perturbation theory), and c is the 
speed of light. 

• Alternatively, by directly matching the couplings of the effective fluid to the measured 
power spectrum, we obtain cl^„^^{A = 1/3) ~ 0.9 x 10~^c^ in remarkable agreement with 
the direct measurement from A^-body simulations. 

• The fluid parameters carry A dependence (as does any 'bare' parameter in an interacting 
field theory). This cutoff dependence is taken to cancel against the cutoff dependence 
of the loop integral. As usual in effective field theories, we 'renormalize' the theory 
by sending the cutoff A — ?> oo and carefully changing the fluid parameters so that 
predictions at low wavcnumbers arc not changed in the process. The finite values of the 
fluid parameters such as c^^^^^ in the A — )■ oo limit is a direct measure of the irreducible 
finite error made in standard approaches that approximate the dark matter on large 
scales as a pressureless ideal fluid. This is an irreducible error that is not recovered 
even by solving non-linearly the equations for a pressureless ideal fluid, as the various 
perturbative approaches attempt to do. This occurs simply because the equations they 
solve are not correct. Our approach, in contrast, should reach arbitrary precision, at 
least in principle, up to non-perturbative corrections. 

• The pressure and viscosity dampen the power spectrum by acting in opposition to 
gravity, which makes sense intuitively. This is able to help explain the observed shape of 
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Ccomb(A = V3) = 0.96 ± 0.1 X 10-6 (c2). 
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the baryon-acoustic-oscillations in the power spectrum relative to standard perturbation 
theory (SPT). 

• More precisely, at one- loop the density-density power spectrum receives a correction 5P 
from the fluid parameters, which we flnd to be 



where Pii,i{k) is the linear power spectrum. Since this is negative and grows as a 
function of /c, the power spectrum is reduced compared to SPT at high /c's, improving 
the agreement with the full non-linear spectrum. 

• We will find that already at one-loop, the computed power spectrum agrees at percent 
level with the non-linear one up to /c ~ 0.24/i Mpc~^ This suggest that in large scale 
structure surveys we should be able to extract primordial information all the way to 
at least such an high wavenumber, improving greatly with respect to the CMB our 
knowledge of the origin of the universe. 

During the years there has been a very large and relevant amount of work in understanding 
perturbatively the large scale clustering of dark matter. An incomplete sample of these works 
is given by [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. 

2 From Dark Matter Particles to Cosmic Fluid 

We take dark matter to be fundamentally described by a set of identical coUisionless classical 
non-relativistic particles interacting only gravitationally. This is a very good approximation 
for all dark matter candidates apart from very light axions. Given that on large scales baryons 
follow dark matter, we can include them in the overall dark matter description. As we discuss 
later, we also neglect general relativistic effects and radiation effects. In this approximation, 
numerical A'"-body simulations exactly solve our UV theory. The coefficients of our effective 
fluid can therefore be extracted directly from the A^-body simulations, following directly the 
procedure described in [4]. Here, the UV theory is described by a Boltzmann equation. 
Therefore, in order to be able to extract the fluid parameters from A^-body simulations, we 
need to derive the fluid equations from the Boltzmann equations and subsequently express 
the parameters of the effective fluid directly in terms of quantities measurable in an A^-body 
simulation. This is the task of this section. 

2.1 Boltzmann Equation 

Let us start from a one-particle phase space density fn{x,p) such that fn{x^'p)d^xd^p repre- 
sents the probability for the particle n to occupy the inflnitesimal phase space volume cPxcPp. 
For a point particle, we have 
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The total phase space density / is defined such that f{x,'p)d^xd^p is the probability that 
there is a particle in the infinitesimal phase space volume d^xd^p: 

n 

We define the mass density p, the momentum density tt* and the kinetic tensor cr*-' as 

pix, t) = -f d'p fix,'p) = -^Yl ^''^("^ - ' (7) 
(X } a 

n 

n\x,t) = -f d'pp'f{x,p) = ^<5(3)(f _ f„) , (8) 
a / a 

a'\x,t) = —J d'ppVfix,p) = -,Vnvi5^'\x - Xn) . 

n 

The particle distribution /„ evolves accordingly to the Boltzmann equation 

Dt at ma^ ox ^ ox op 

where 0„ is the single-particle Newtonian potential. There are two important points to 
highlight about the former equation. First, we have taken the Newtonian limit of the full 
general rclativistic Boltzmann equation. This is an approximation we make for simplicity. 
All our results can be trivially extended to include general rclativistic effects. However, it 
is easy to realize that the Newtonian approximation is particularly well justified. Non-linear 
corrections to the evolution of the dark matter evolution are concentrated at short scales, 
with corrections that scale as k'^ /k%j^. General rclativistic corrections are expected to scale 
as k^/{aHy. This means that we should be able to cover up to wavelength of order 300 Mpc 
before worrying about per mille General relativity corrections. Furthermore, one of the main 
goals of this paper is to recover the parameters of the effective ffuid of the universe from very 
short scale simulations valid on distances of order of the non-linear scale. The parameters we 
will extract in the Newtonian approximation are automatically valid also for the description 
of an effective fiuid coupled to gravity in the full general rclativistic setting. 

A second important point to highlight in the former Boltzmann equation is about the 
single-particle Newtonian potential 0^. Following [4], the Newtonian potential is defined 
through the Poisson equation 

d'^(f)^4TiGa^p- Pb) , (10) 

with pb being the background density and 9^ = S^Widj. We raise and lower spatial indexes 
with Sii. The solution reads 



(p = 2_^(pn+ ^ , (11) 



p- 

Gm 



{x) = - e-^l^-^"l . (12) 
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Notice that the overall is IR divergent in an infinite universe. This is due to a breaking 
of the Newtonian approximation. We have regulated it with an IR cutoff /i that we will take 
to zero at the end of the calculation. Our results do not depend on /i, as indeed we are 

interested in very short distance physics. 

By summing over n, we obtain the Boltzmann equation for / 

Rl = ^ + J_.^_rn y ^.^ = (13) 
Dt dt ma? dx dx dp 

n,n;n^n 

2.2 Smoothing 

Following [4] , we construct the equations of motion for the effective fiuid by smoothing the 
Boltzmann equations and by taking moments of the resulting long-distance Boltzmann equa- 
tion. The smoothing guarantees that the Boltzmann hierarchy can be truncated, leaving us 
with an effective fluid, indeed, notice that it is not trivial at all that we should end up with an 
effective fluid. Fluid equations are usually valid over distances longer than the mean free path 
of the particles. But here for dark matter particles the mean free path is virtually inflnite. 
What saves us is that the dark matter particles have had a finite amount of proper time, 
of order H~^, to travel since reheating, and they traveled at a very non-relativistic speed. 
This defines a length scale vH~^ ~ l/kNi which is indeed of order of the non-linear scale. 
This length scale plays the role of a mean free path, as verified in [4] . The truncation of the 
Boltzmann hierarchy is regulated by powers /c/Zcjvl ^ 1- 
We define the Gaussian smoothing 

W^{x) = ("^Ve-i^'-' , W^{k) = e-^x^ , (14) 



with representing a /c-space, comoving cutoff scale. This will smooth out quantities with 
wavenumber A; > A, or equivalently with waveleghts smaller than A < 1/A. We regularize our 
observable quantities 0{x, t), p, tt, 0, . . . , by taking convolutions in real space with the filter, 
defining long-wavelength quantities as 

Oi{x,t) = [0]^{x,t) = J d^x'WA{x-x')0{x') . (15) 

Notice that in Fourier space W{k) — > 1 as A; — > 0: our fields are asymptotically untouched at 
long distances. 

The smoothed Boltzmann equation becomes 



Dt 



'^'-£.-i-"E./^VH.4.--.^)^(.^).f . (16, 



. dt ma^ dx _ 

■'^ n,n,n^n 



Fluid equations are obtained by taking successive moments 

'Df 



drp p'^ . . . p*" 



Dt 



(i;,p) = 0, (17) 

A 
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creating in this way a set of coupled differential equations known as Boltzmann hierarchy. 
As we will explain in more detail later, it will be sufficient for the purposes this paper to 
stop at the first two moments (one-loop approximation). The first two moments will give the 
continuity and momentum equations in the approximation in which the fiuid is described by 
the Navier-Stokes approximation, with the addition of a stochastic term. We obtain 

pi + mpi + -di{piv\) = , (18) 
a 

v\ + Ev\ + -v\d,^, + -d,<^^ = --^9, [r'^]^ . (19) 

yjj (Jj CL tJl 

Let us define the various quantities that enter in these equations. We define the long wave- 
length velocity field as the ratio of the momentum and the density 

= ^ ■ (20) 
Pi 

The right hand side of the momentum equation (19) contains the divergence of an effective 
stress tensor which is induced by the short wavelength fluctuations. This is given by 

\r'']^ = x^ + ^ , (21) 
where k and $ correspond to 'kinetically-induced' and 'gravitationally-induced' parts: 

Kf = - piv\vi , (22) 
1 

■ 871^02 



where 



(23) 



Note that we have subtracted out the self term from w]-^ , as necessary when passing from 
the continuous to the discrete description in the Newtonian approximation, and used that 
d^(j) = 47rG'a2(p — p^) and (9^0; = AnGa^i^pi — pb) to express in terms of (p and (f)i. In the 
limit in which there are no short wavelength fluctuations, and A oo, ki and vanish. 
In App. A we provide the above expression written just in terms of the short wavelength 
fluctuations. 



2.3 Integrating out UV Physics 

The effective stress tensor that we have identified is explicitly dependent on the short wave- 
length fluctuations. These are very large, strongly coupled, and therefore impossible to treat 
within the effective theory. When we compute correlation functions of long wavelength fluctu- 
ations, we are taking expectation values. Since short wavelength fluctuations are not observed 
directly, we can take the expectation value over their values. This is the classical field the- 
ory analog of the operation of 'integrating out' the UV degrees of freedom in quantum field 
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theory, now applied to classical field theory. The long wavelength perturbations will affect 
the result of the expectation value of the short modes, through, e.g., tidal like effects. This 
means that the expectation value will depend on the long modes. In practice, wc take the 
expectation value on a long wavelength background. The resulting function depends only on 
long wavelength fluctuations as degrees of freedom. In this way, we have defined an effective 
theory that contains only long wavelength fluctuations. Since long wavelength fluctuations 
are perturbatively small, we can Taylor expand in the size of the long wavelength fluctuations. 
Schematically we have 

Si + ... . (24) 



For the precision wc pursue in the rest of the paper, we will stop at linear level in the 
long wavelength fluctuations, though nothing stops us from going to higher order. By the 
symmetries of the problem, the resulting stress tensor must take the following form 

(25) 

This is the stress tensor of an imperfect fluid, pb is the background pressure that is induced 
by short distance inhomogeneities even in the absence of long wavelength fluctuations, 
is the speed of sounds of the fluctuations: Sp — c^Sp. The parameters and Csv are the 
coefficients for the bulk ( and the shear r) viscosity respectively, with units of velocity. They 
are related to rj and ( by the relation r] = 3pfeC^^/(4i7), ^ = pbC^^/H . Ar'-' represents a 
stochastic term, that takes into account the difference between the actual value of r*-^ in a 
given realization and its expectation value ^. We will come back to this term shortly, but 
it is worth noting that neglecting this term in the above equations reproduces the familiar 
Navier-Stokes equations. 

Finally, the ellipses (...) represent terms that are either higher order in 5i, or higher 
order on derivatives of 5i. Indeed, higher derivative terms will be in general suppressed by 
k/kisfL -C 1, and, as typical in effective field theories, wc take a derivative expansion in those. 
Astrophysically, these terms would corresponds to the effects induced by a sort of higher- 
derivative tidal tensor. Once we expand in derivatives of the long wavelength ffuctuations, 
we take the parameters in (25) to be spatially independent, but time dependent. 

The coefficient Spb, Cg, Csb, Csv are determined by the UV physics and by our smoothing 
cutoff A, and are not predictable within the effective theory. They must be measured from 
either A'"-body simulations, or fit directly to observations. This is akin to what happens 
in the Chiral Lagrangian for parameters that can be measured in experiments or in lattice 
simulations, such as F^^. We first define the correlation functions that will allow us to extract 
these parameters from small A^-body simulations. 

^For the readers familiar with the in- in formaUsm, this term wiU take into account the cut-in-the-middle 
one-loop diagrams [26]. 
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2.4 Matching Correlation Functions 

It useful to define the following quantities from the stress tensor 

Ji^-^d,[T^^]^ , Af^U'Jl, (26) 

(J/L/fj (Jj 

and to introduce a dimensionless velocity divergence 

Then, according to (25), we have 

aj; = c'ASi + IcldjQf + (t + ^l) ' (28) 

a^A, = cld% + {cl + cl) d-Q, , 
a'B, = cld'Q, . 



In order to extract the parameters of the effective fluid, we multiply each of these functions 
with long wavelength fields, and take expectation values. By forming suitable combinations 
of these, the parameters of the effective fluid can be extracted. We will need the following set 
of correlation functions 

Pa5{x) = {Ai{x' + x)Si{x')) , (29) 

PAsix) = {Ai{x' + x)Qi{x')) , 

Pa^^bJx) = {Al^x' + x)ei,,ix')) , 

Pseix) = {Biix' + x)Qi{x')) , 

Pss{x) = {6i{x + x)6i{x )) , 

Pseix) = {6i{x' + x)&i{x')) , 

Peeix) = {Qiix' + x)eiix')) , 

Pe.«0^(^) = (ef(f' + f)ef,(a?')) , 



where Si = Spi/ph. From these we obtain the following expressions for the parameters of the 
effective theory 

^2 ^ ^2 PAe{x)d^P5e{x) - PA5ix)d^Pee{x) 
{d^Pse{x)f -d^P56{x)d^Pee{x) ' 

2 ^ 2 PA5{.x)d'^P5e{.x) - PAe{x)d'^Psi{x) 
''''''' {d^Pse{x)f -d^Ps6{x)d^Pee{x) ' 

2 ^ 4 2 PA^'-euii.^) - PAe{x) ^ 2 PBejx) 
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where = c^^ + c^^ is the sum of the viscosity coefficients. By extracting the correlation 
functions in (29) from A'"-body simulations, and performing the ratios in (30), we should be 
able to extract the parameters of the effective theory. Notice that the ratios are supposed 
to be spatially independent. Such behavior is expected to hold at large distances x ^ 
where the higher derivative terms are negligible. 

In simulations, we should in principle also measure the stochastic components of the stress 
tensor. In the two point function at one-loop, at leading order in derivatives, it enters just 
the correlation function of the trace. This amounts to measuring 

(j;(f,t)j/(f+f,0) ■ (31) 

We will see that the effect of this stochastic term is accidentally higher order in 5i and so does 
not enter at leading order. 

After all these parameters have been measured in A^-body simulations, the EFT is prone 
for perturbation theory. It is alternatively possible to perform directly perturbation theory 
and fit the results to observations. We will be able to perform both approach and check that 
we obtain the same result. We will describe in detail how to measure these quantities in 
simulations in App. D and we will give the results of these measurements in Sec. 4. For the 
moment we will instead directly move to apply perturbation theory with our EFT. 

These parameters can either be measured from A^-body simulations directly or kept generic 
and then extracted by fitting the results to observables. We are able to do both and verify 
that we obtain the same result. We will describe in detail how to measure these quantities 
in simulations in App. D and we will give the results of these measurements in Sec. 4. First 
we develop perturbation theory within the EFT that will allow us to make predictions and 
to extract these parameters from observations. 



3 Perturbation Theory with the EFT 

We now proceed to perform perturbation theory within our EFT. The non-linear equations 
of motion that we need to solve are 

V'ct>i = lHln^''\ + ... , (32) 

5i^--d,{{l + 5i)vf) , 

Qj 

where ' = d/dt, H = d/a, A J* = dj{AT^')/{apb), flm is the present day matter fraction, 
oo is the present day scale factor, usually taken to be equal to 1, and . . . represent higher 
order terms (in sense that we will explain shortly) that come from the expression of the 
short wavelength stress tensor r'-' in terms of long wavelength fiuctuations. Our theory is 
defined on scales longer than the non-linear scale. For this reason we have 5^ <^ 1. The long- 
wavelength velocity vi and the long-wavelength (pi are small even inside the the non-linear 
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scale and they are even smaller at larger distances. This means that we can reliably solve the 
above non-linear equations iteratively around the linear solution. Such an iterative solution 
is very similar to what is done in quantum field theory, where the solution to the quantum 
non-linear equations is organized in Feynman diagrams. Indeed we can organize the various 
perturbative terms around Feynman diagrams even in this case. The result is very similar to 
what is computed in the in-in formalism, for example when computing quantum corrections 
to inflationary correlation functions [26]. Indeed, the calculation we are going to do shares 
many of the features that are present in normal quantum field theory computations: cutoff, 
renormalization, running, and so on are all concepts that will appear and prove useful as we 
proceed. They have nothing to do with the word 'quantum' in 'quantum field theory', rather 
they have to do with the 'field theory'. Our calculation is for a classical field theory and 
shares all these features. 



3.1 Organization of the perturbation theory 

The simplest way to organize our perturbation theory is to use the fact that, in any order of 
magnitude approximation, is constant at all scales, of order 10~^, and that well inside the 
horizon the Newtonian approximation holds. For length scales longer than the equality scale, 
at the linear level we therefore have 

0, (Ax ~ L) ~ 10-^ , vi{Ax L) r^lQ-^^ , (33) 

Si {Ax r^L)r^ 77^0/ (Ax - L) - 10" " ^ 



^2^2-^' V / ^2^2 

We see that as L ^ 0, 5; grows and indeed becomes of order one at L ~ Aat^ ^. At distances 
larger than the non- linear scale, we therefore expand in powers of 5;, keeping in mind that 
the additional fiuctuations scale as in (33). Let us estimate the relative size of the terms. 

Loop corrections: it is easy to estimate from the non linear structure of the equations 
that diV^ ~ H5i. Notice that in the power spectrum we need to take two non-linear corrections, 
or alternatively look at the cubic corrections. The non-linear terms scale as 

non — linear terms SiV'ldjV] H Sfv] 

~ — <^^L± M- ~ Sf (34) 

Hubble friction Hv} Hv\ ' ^ ' 

Loop corrections therefore scale as 5f , peaked at the highest possible scale within the theory A. 

Pressure and Viscosity terms: these terms result from integrating out the modes 
higher than the A scale. So, naively they should scale as the 5 at the extreme UV scales 
beyond the effective theory. Since the theory in the UV is strongly coupled, very large 



^For the propose of estimating at order of magnitude level, we have assumed that the k modes under 
consideration are longer than the equality scale fccq ^ 0.01 Mpc^^ and taken (6; to be fc-indepcndent. On 
shorter scales (j)i decays, and the estimates need to be slightly modified. This subtlety will be important 
for the actual numerical contribution of the various terms, and it will be properly accounted for, but is not 
particularly relevant for the order of magnitude estimates that control our power counting, and so we will 
ignore it for simplicity's sake. 
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corrections are expected, and the result cannot be extrapolated from the linear regime, even 
at the order of magnitude level. This is why we will measure the parameters such as from 
A'^-body simulations. What wc will find is that these parameters are of order 10~^. This 
happens because the combination of short modes that generates the parameters like Cg are 
such that short modes that have virialized do not contribute ^. Since these terms scale as (f)6, 
the contribution is peaked at those modes that have just become non-linear 5 ~ 1, but not 
yet viriahzed. We therefore expect the parameters c^'s to be of order ~ 10~^. We will see 
that the fact that short modes entered the horizon in the radiation era makes this a bit of an 
overestimate. At this point we are ready to estimate the size of these corrections: 

Pressure , Viscosity ctddi o^'^^i , . 

~ ~ c ~ oi . (35) 

Hubble friction Hv} ' H^i IQ-^ ' ^ ' 

Notice that thanks to the strongly coupled UV theory (or thanks to the fact that non-linear 
structures virialize), we have that for c\ ~ 10"^, the contribution from these terms is larger 
than the one loop contribution in the low energy theory. This is so because the theory is 
strongly coupled in the UV. We conclude that one insertion of these terms counts at least as 
a one-loop term. 

Stochastic terms: Let us now continue on to evaluate the effect of the stochastic terms. 
This is a bit more complicated. Let us evaluate the relative effect on the power spectrum. 
The structure of the equations leads to the following approximate non-linear solution 

2 d'^ (9^ Ar 

Sl,non-lin. ~ Sl,lin + Cgj^Si^Un + ■ (36) 

In the power spectrum we therefore have 

(<^.<^.)i-ioop - cl^,{Sf) + ' {At') , (37) 

as the stochastic part must be correlated with itself. Due to virialization, we expect that the 
correlation function of r should be Poisson like on independent pixels of order the non-linear 
scale k]^\. We therefore estimate 



<Ar^) ~ ( j*^) ~(c?p.)^f^) . (38) 



kNL J \kNL 

This is indeed confirmed by calculations in perturbation theory [4] . We therefore have 

Stochastic k Stochastic A; , ,s/2 / x 

-Si - 5f/' . (39) 



Pressure k^L ' Friction k^L 



This tells us that on scales longer than the non-hnear scale, the contribution of the stochastic 
pressure is parametrically smaller than the pressure effects. Since in this paper we will stop 



^This was found and used in [4] to show that there is very little backreaction on the evolution of the 
universe from short scale non-linearities. Short scale gravitational collapse changes the equation of state of 
the overall universe by a relative factor of order 10~^. 
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at one-loop order, we can therefore neglect this correction. It should be noted that these con- 
tributions scale parametrically differently than the loop contributions, and so, depending on 
the scale considered, they might be more relevant that a 2-loop contribution. We notice that 
something similar happens also at the level of dissipativc fluids, where we generically include 
dissipativc terms through the Navier-Stokcs equations, but we neglect stochastic terms. 

Higher Derivative Terms: When we take the expectation value of the short-distance 
stress tensor in a background of a long mode, we have Taylor expanded both in the long 
wavelength fluctuations and in their derivatives. Higher power corrections scale with powers 
of 5i, and so corresponds to higher- loop terms. Higher derivative terms instead scale as powers 
of k'^/kffL ^ where we have taken the squared because of rotational invariance. This shows 
that higher derivative terms scale nicely as loop terms. If we allow the cutoff to remain finite, 
we should also include higher derivative terms that scale as fc^/A^. 

General Relativistic and Radiation Corrections: In this paper we will neglect gen- 
eral relativistic corrections and all non-linear contributions coming from the fact that the 
universe was radiation dominated at early times. General Relativistic corrections scale as 



For the high scales where non- linear corrections are relevant, for example at the BAO scale 
kBAO ~ 10~^, these corrections are of order 10~^, and so uninteresting from this point of view. 
Radiation is the dominant component of the universe at early times. Neglecting it amounts 
to neglect corrections that scale as a/agq ~ 10~^, where a is the scale factor and agq is the 
scale factor at matter radiation equality. Inclusion of these corrections in perturbation theory 
has been studied in [20], and it gives a small correction to power spectrum, and small, but 
potentially measurable, corrections to the three-point function, corresponding to /jvl ~few. 
Both General Relativistic effects and radiation effects do not represent an intrinsic limitation 
of our EFT. They can be straightforwardly included in our formalism, by simply improving 
the equations of motion we use in this paper. 

In summary, we see that apart from the stochastic terms and some higher derivative 
terms, all the remaining terms: loops, pressure, higher derivatives and higher powers of 6i 
from T, scale as powers of 5i, which is our main ordering parameter. Stochastic terms instead 
contribute at leading order as 5f^^, so they count as one loop and a half. Cutoff-dependent 
higher derivative terms scale as /c^/A^. 

Cutoff dependence and effective expansion parameter: So far, it looks like that 
our expansion parameter is the highest Si we have in our theory, which is 6i{k ~ A) ~ 
A'^/k%]^. However, the situation is even better than this. So far, we have defined our theory 
with a regulating cutoff at A; ~ A. Because of this, all our intermediate results depend 
exphcitly on A: Cs{A), Cs^(A), etc. and loops need to be cutoff at A. This induces an exphcit 
A dependence plus higher derivative terms of order /c^/A^. However, the sum of all the 
diagrams will be independent of A. Indeed Cs(A), etc. should be really thought of as one- 
loop counterterms. Upon carefully choosing the counterterms Cs(A), etc., any A dependence 
cancels apart from terms in A;/ A that should be removed by higher derivative corrections in 




(40) 
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the stress tensor that we neglected. In order to resolve such error with the least effort, we 
will choose the counterterms at a fixed cutoff in such a way as to have the theory agree with 
observations at a certain renormalization scale /cren.j and then we will extrapolate our results 
to A — 7- oo, effectively letting the residual terms in /c/A vanish. In this A — )■ oo regime, the 
loop term gets dominated by the regime in which one of the modes has wavcnumber of order 
the non-linear scale, while the other has a wavenumber of order of the external wavenumber. 
In this way, loop terms scale as one-power of 5i as the counterterm ^. At this point, the 
expansion parameter of the EFT will be Si ~ k^/kj^^ evaluated at the scale of the external 
modes, with no residual A dependence even in the expansion parameters. 

Again, this is very similar to what happens when one computes loop corrections in the 
Chiral Lagrangian. After regulating the chiral theory with a cutoff A, there are naively two 
expansion parameters. If E is the energy scale of the process, we have E / F^^ and i?/A. After 
renormalization and by sending A ^ oo, we are left only we E/ F^^ as an expansion parameters. 

In summary, the expansion parameter of the EFT is ~ (k/kNL), where k is the typical 
wavenumber of the external modes. Loops in the EFT, counterterms and higher-derivative 
terms scale as Si. Stochastic terms start contributing at order 5/ . 



3.2 One- loop Perturbation Theory 

We are now ready to implement perturbation theory for the power spectrum at quartic order 
in Si, that is at one-loop. At this order, the equations we are going to solve are the ones 
in (32) with A J and the . . . terms neglected. 

Let us write the equation for the vorticity wl = e^^^djVk- Neglecting the stochastic terms 
that we argued are small, we have 

In linear perturbation theory the vorticity is driven to zero, and this occurs even the more 
so at this order in perturbation theory, as the source is proportional to Wi. While at higher 
order one could expect vorticity to be generated, at this order, and therefore for the purposes 
of this paper, we can take it to be zero. This means that we can work directly with the 
divergence of the velocity 

Oi = div\ (42) 
Using a as our time variable, the equations (32) reduce to 

a-HS'i + ^o^io'^ k - q)Si{k - q)ei{q) , (43) 

ane[ + nei + I'^Si - clk\ + '^Oi^-J -^,m k - mfi - mio) , 

^Another kinematically allowed possibility is for the modes to have both wave numbers close to the non- 
linear scale, but slightly different so that their sum is equal to the external wavenumber k. This contribution 
would naively scale as ((5;)° ^ 1. However this contribution in this regime scales as the stochastic term df^^. 
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where l-L — a ^da/dr, subscript o for a quantity means that the quantity is evaluated at 
present time, we have set oq = 1, ' represents d/da and 

(k + q] -k (k + q\k-q 

a{k, q) = ^ ■/ , m q) = ^—{f-^ ■ (44) 

As we discussed, the parameters Cg, and Csv are time dependent and must be measured 
in the simulations as a function of time. For the purposes of this paper, we will make the 
simplifying assumption that their time dependence can be inferred in perturbation theory. In 
other words, we will measure them at one time and deduce their values at different times by 
perturbation theory ^. 



3.2.1 Perturbative Solutions 



Since the correlation function of matter overdensities is small at large distances, we can solve 
the above set of equations (43) perturbatively in the amplitude of the fluctuations. For 
the computation of the power spectrum at one loop, it is enough to solve these equations 
iteratively up to cubic order. Order by order, the solution is given by convolving the retarded 
Green's function associated to the linear differential operator with the non-linear source term 
evaluated on lower order solutions. At second order we obtain 

ddG{a,d)d^n^{d)D'{d)A [2 [ d^q/3{q,k - q)Ssi{k - q)5si{q) 



+ (^J^ddG{a,d) (^2d^H\d)D'{df + 3nln^^^^ 



X 



d'^qa{q, k — q)5si{k — q)Ssi{q) 



Let us explain some of the relevant expressions that appear here. G{a, d) is the retarded 
Green's function for the second order linear differential operator associated with S that is 
obtained after substituting 6 in the second equation of (43) with the value obtained from the 
first, and linearizing. In doing this, it is important to neglect all the terms of order because, 
in our power counting, they count as non- linear terms. This is given by 

-a^n^{a)dlG{a, d) - a {2H\a) + aH{a)H'{a)) daG{a, d) + 3^^^G(a, d) = 5{a - d) , 

G{a, a) = for a <d . (46) 

^As we will see, these parameters need to cancel the A dependence associated to the regularized loops. 
The part of these parameters that depends on A can be therefore reliably inferred in perturbation theory. 
However, the part of these parameters that is A-independent and that represents the finite contributions 
should be measured in simulation or in observation. We will assume that the time dependence for these two 
components is the same. We will check that this is an accurate approximation in an upcoming paper [27]. For 
this approximation, we stress that since these are 1-loop terms, it is important to know them up to a relative 
factor of order 6i <^1. 
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For a ACDM cosmology the result can be expressed^ as a hypergeometric function, although 
its form is not particularly illuminating. For all calculations presented here it is sufficient 

to numerically solve the above differential equation. This can be easily accomplished by 
replacing the d{a — a) on the RHS of the first equation with zero, but starting with the 
boundary conditions being G{a,d)\a=a = 0, and ^G{a,a)\a=a = l/(a'H(a))^ . fn principle, 
it is possible to include in the linear equations that determine the Green's function and the 
growth functions also the higher-order linear terms proportional to and c^. Doing this 
amounts to resumming the effect of these pressure and viscous terms. The resulting linear 
equation can be easely solved numerically, finding for example that the growth factor becomes 
/c-dependent, being the more suppressed the higher is the wavenumber [21]. However, it is 
not fully consistent to resum these terms without including the relevant loop corrections. 

D{a) represents the growth factor at scale-factor-time a. In particular, we have written 
the linear solution as 

4'\k,a) = ^6s,ik), (47) 

with qq being the present time, and Ssi representing a classical stochastic variable with 
variance equal to the present smoothed power spectrum 

{5s^(k)5s,(q)) = (27r)'S^'\k + ^PnAk,^) , (48) 
with Pii^i{k) being the smoothing of the linearly computed power spectrum at present time 

Pn,i{k,^) = WAikyPnMk) ■ (49) 

A very useful simplification is due to the fact the growth factor and the Green's function 
are /c-independent. This is due to the fact that at linear level we can neglect the pressure and 
viscosity terms that would otherwise induce a fc-dependence. Because of this, the convolution 
integrals that would couple time integration and momentum integration nicely split into 
separate time integrals and momentum integrals that can be simply performed separately. 
We have tried to underline this in (45) by adding suitable parenthesis. Iterating, we obtain 
the solution for S at cubic order 6^^\ For brevity, we report it in App. B. Notice that in the 
terms in S^^^ (and 5^^') we have neglected the contribution from the pressure and the viscosity, 
which count as third order terms. They give 

-> k'^ f"" -> 

^i%omS^' ^ ~D{^ Jq *^^"' D{a) 6si{k) , (50) 

where cl^^-^ is given by 

Ccomb(a) = c,{a) + «^^^c^(«) , (51) 

and it is the combination that is relevant at one-loop order. For the terms multiplying and 
cl in the second equation of (43), we can substitute the linear relation 

9l^\a, k) = -anda5f\a, k) = -an^^6l'\a, k) . (52) 



^ Using, e.g., Mathematica's "DSolve" function. 



17 



Notice that 5^2) and S^^^ are the same as the standard ones used in perturbation theory, 
with just two differences. The first is the important smoothing of the sourcing power spectrum. 
This makes the convolution integral that wc arc going to perform next rapidly converging, 
but also A dependent. The second difTerence is of a more technical nature, and relies on the 
fact that in standard perturbation theory the time dependence of the non-linear solution is 
approximated by the linear growth factor D elevated to the power 2 and 3 for 6^"^^ and d^^\ 
while the momentum dependence is approximated to be the same momentum dependence as 
in standard EdS universe. This procedure is exact in EdS, but not so in other space times. 
Some studies [18] (see also [29, 30, 31]) have checked that this is correct up to percent level 
on the full power spectrum. Since however percent accuracy is the target of next generation 
experiments, we decide to perform the correct computation, which is not so very complicated 
to set up in any case. For the purpose of comparing with the literature and to gain familiarity 
with the EFT setup with simpler formulas, we provide results obtained with this approximate 
treatment of the perturbed solutions in App. C. 



3.2.2 Diagrams 

By contracting the non-linear expression we obtain the non-linear corrections. There are three 
diagrams at order 5f. After including the linear contribution, we have 

{Si{k,ao)Si{q,ao)) = {27r)^ S^^\k + q) (Pn(A;, oq) + P22(A:, oq) + Pi3(A:, ao) + P13, ,(A;,ao)) 

(53) 

with 

Pn{k,ao) = {S^'\k,ao)S^'\q,ao)y , (54) 

P22ik,ao) = {6l'\k,ao)5l''\q,ao)y , 

Pi3(A;,oo) = 2{6i''\k,ao)5i'\q,ao)y , 

As, ,(^, ao) = 2(<5f^l {k, ao)Sl'\q, ao))' , 

' comt. ''"-comb 

where the (. . .)' means that we have removed a factor of {27r)^S'^^\k + q) from the expectation 
value. Pii represents the unsmoothed linear power spectrum, as the linear theory does not need 
to be regularized. The term P13 ^2 ^ is supposed to remove the A dependence that comes 
from P13. It is a countcrterm diagram. Strictly speaking, we would need a countcrterm 
diagram also from P22, which is provided by the two-point function of the stochastic source 
A J* in (32). As we discussed, this term is supposed to count as a 5f term, and therefore we 
neglect it. This means that the A dependence associated with P22 is very weak at this order 
in the calculation. The full stochastic term will be included in a following paper [27]. 
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The expressions for P22, P13, Pis are given by 

' comb 



P22ik,ao) 



dq d(cos 9) - 



,Pn,i{Q,^)Pii,i{\k-q\,A-) 



ao 



167r2D(ao)4 J " ' {k^ - 2 cos{9)k q + q'f 

da a?G{aQ,d) n'^{d)D'{df^ Acos{0){k - cos{0)q) 

+3nlnm U da G{ao, a) j (cos(^) {k - 2 cos{9)q) + q) 

da d'^G{ao, a) 'H\d)D'{d)A (SA;^ cos(^) - A; g (4 cos(^)2 + l) + 2 cos{9)q'^) 



m 

where 008(6*) = k ■ q/{kq); 

Pi3{k,ao) = 



ao 



da G{ao, a) ) (^^ - 2 cos(^)A; q + q^) 



(55) 



2 e 



96{27ryD{ao)^ 
dr 

/ — [l2r^V4 - 24rP5 + 4r=' (I6P1 + 8P2 + 4^3 - 3^4 + 24P5) 

+8r^ (4D2 + 2V3 - 6D4 + 3D5 - 4I?6) 



(56) 



+3 {r^-lf{r^V^ + 2Vr,) log ( ^ 



l-r) 



(1 + r)^ 



where 2^i,...,6 are given in App. B, and finally 

Jf2 /•ao 



D{ao) Jo 



da G{ao, a) c^„^^{a) D{d) Pn,i{k, A) . 



(57) 



The convolution integrals in P22 and P13 are the sign that these are one loop diagrams. 
P13 „2 does not have a convolution integral as it is a one-loop counterterm. The diagrams 

' comb 

— * — * — * — * 

are pictorially represented in Fig. 1. Since a{k, —k) — (5{k, —k) — 0, there are no non-lPI 
diagrams. Notice that in P13 and P22 we have already carried out at least a part of the 
angular integration. The sign of the P10 /2 in (57) is also quite intuitive. For positive Cg or 

' ' comb 

Cy, the gravitational collapse is slowed down, and so this contribution tends to decrease the 
gravitational collapse. 



3.2.3 Cutoff- (in) dependence 

Each diagram is dependent on the cutoff A: P22 and P13 through the smoothed linear power 
spectrum, while Pi3^ ^^ depends on A through c^comh^ which is A dependent because it 
arises from integrating out the short distance fluctuations. The A dependence of P13 is to be 
cancelled by P13 ^2 , while the one of P22 from the stochastic fluctuations. 

' comb 

Consider the sum of the P13 and P13 terms. In order for this sum to be A independent, 

' comb 

both must have the same A;-dependence in the relevant regime. By inspection of (57), we see 
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Figure 1: Diagrammatic representation of P22 (top left), P13 (top right), and Pio 2 (bottom). 

' comb 

Continuous green lines represent Green's functions, red dashed lines represent free fields, and red 
crosses circled by a dotted blue line represent correlation among free fields. 



that c2 goes as fc^, which imphes that P13 should behave in the same way in this 

' comb 

regime. This is in fact the case, as can be readily verified by taking the k limit of (56). 
In particular, we can define a A-independent renormalized parameter Ccomb;ren. defined at a 
renormalization scale k^en. and a A-dependent counterterm parameter c^omb;ctr. (^) ^^ch that 

'^comb('^) ^) ~ '^comb;ren.('^) ^rcn.) + C(,ojjj]-,.f,tj, (a. A) . (58) 

Ccomb ctr.('^' ^) Kiust havc the same time and A-dependence of P13, while Ccomb;ren.(a, ^ren.) 
is determined by matching to simulations or to observations at a specific k = /Cren.- The 
time dependence of c^omb;ren.('^) could in general be different from the one of c^omb;ctr.(^' 
However, for the purposes of this paper, we can approximate them to be equal, and we will 
check this approximation in a forthcoming paper [27]. We can therefore extract the time 
dependence of c^Qj^b(a, A) from the A; — )■ hmit of P13. We obtain 

2 2 -^^2 (a) 

'^comb;ctr.('^' ^) = combictr. ('^0; ^) ^ ''°™''/ T ) (59) 

comb 
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where 



D^2 (a) 



[{GSn^oQlTiia) + UnlQmMa) + h2a^'H^{a)D' {af) (60) 



D'ia 



and 



7i(a) 



da G{a, a) — z — , 



da daG{a, a] 



a 

D{af 
a 



FA a) 



f ddG{a,d)d^D'{dfH^{d), (61) 
Jo 

da daG{a,d) d^ D\dY'H'^ [a) . 



A plot of the time dependence of the speed of sound is given in Fig. 2. 

Time evolution of c^comb (A=oo, kren= ■ 16) 




Figure 2: Time dependence of c^^^^^i^ as inferred using the correct time dependence from P13 and 
instead using the approximate time dependence derived from the growth functions (see App. C). 
Starting from very early times, we see that c^^^^^ grows as a functions of time, peaks at about 
a ~ 0.7, and then decreases near the present epoch, probably as due to the onset of dark energy. 



'comb 



is positive, implying that this term tends to slow down the collapse of structures. 



We determine the value of c^ombxtr.('^O) A) in two different independent ways: one involving 
fitting to observation of the power-spectrum derived from simulations, and the other involving 
direct measurement from simulation. In the first, simplest, way, we determine c^oj^b-ctr.('^O) A = 
00) by matching the one-loop EFT power spectrum at k = fcren. to the power-spectrum 
extracted from simulations (or directly from precise observations in the future!). We can 
do this at various values of A, but we take the A — )■ cxd limit in order to drive to zero any 
effect from higher derivative terms down by powers of k/A. At this point, we can derive 
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'^comb('^O) A 7^ oo) by running the value at A = oo down to a finite A. The formula is given by 



Ccomb(«0,A 7^ Oo) 




(1.2 rao 





(62) 



The limit fccxt — >■ is necessary in order to suppress higher derivative terms down by 
powers of fccxt/A. The running of c^^j^b is plotted in Fig. 3. c^^^^^{ao,A = oo) ^ 6.2 x 10"'' 
and is the value obtained by fitting to data using kren. = 0.16 /i Mpc~^ ^. We see that as A 
decreases, we integrate out more and more modes, and 0^^^^^ grows. We see that the result 
matches with the one obtained by the second method for measuring c^comb^ ^^^^ using 
A'^-body simulation to extract directly the correlations in (30). We perform measurements at 
A = 1/3 /i Mpc~^ and A = 1/6 h Mpc"-*^, and we see that the match is extremely good. We 
describe more precisely how these measurements are derived in Sec. 4. We take this as an 
extremely promising indication of the strength of our approach. 

In order to elucidate the effect of the higher derivative terms, we plot in Fig. 4 the value 
of cl^^rahi^o, A = 1/3), for various values of the external /Cext- It is only for vary low /Cext's that 
•^comb becomes /cext. independent. 

Finally, there is a third method in which we could have derived c1^^^{ao, A). By keeping 
A finite, we could have fit the analytical results to iV-body simulations by including higher 
derivative terms proportional to powers of k/A. Indeed, unless A — )■ oo, the largest of these 
terms arc not negligible and need to be included to get the correct c^^^^^^laQjA). A description 
of this approach in detail is given in App. E, and leads to the same results for CpQjjjb(ao, A). 



If this language of effective field theory is to be born out, we must be able to take the 
fundamental theory, integrate away UV effects, and find agreement in terms of the parameters 
described above. Fortunately we have, in the form of simulation, exactly those calculations in 
the fundamental theory. By smearing the positions of simulated particles with a normalized 
Gaussian function of width A, we are able to introduce a soft UV cutoff of order 1/A. For 
correlations on scales longer then the cutoff we can directly measure c^^^^ corroborating the 
perturbative analysis presented above. 

Specifically we consider random downsamples of positions and velocities of the Consuleo 
simulation^ [28] from 2.7 x 10^ particles downsampled to 1, 000, 053 particles distributed over 
(420 Mpc)^. Even this incredibly coarse resolution allows us to measure the following fields 

''It is somewhat interesting to notice that by using Cs„ ~ 10"'' we find a shear viscosity of order ~ 20 Pa s, 
in SI units. This value is very similar to that of some everyday items such as chocolate syrup! 

*The simulation parameters are: Vim = 0.25, JIa — 0.75, h = 0.7 {H = 70 km/s/Mpc), erg = 0.8, and 
Hs = 1 with measurements described taking place at z = 0. 



4 Fluid parameters from A^-body simulations 
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Fi gure 3: Running of c^^^^^ as a function of A. The purple band contains the region for the 
values of c^^^^^ as inferred from matching with the non-linear power spectrum from CAMB at the 
renormalization scale k = O.lh Mpc^^ and k = 0.18h Mpc^^ . The dependence on the renormalization 
scale is a measure of the importance of higher loops. We see that as A ^ oo, c^comh decreases as 
more and more modes are included within the regime of validity of the EFT. However, the fact that 
as A = oo, c^ojjjb 7^ is an indication of the fact that the fundamental theory is not described by 
a pressureless ideal fluid, but by indeed freely streaming dark matter particles. Data points with 
1 cr error bars represent the value obtained from N-body numerical simulations using the methods 
described in Sec. 4 using two different smoothing lengths A^^ . Given the error bars from numerical 
evaluation, the measured values are in remarkable agreement with what inferred from renormalizing 
using the power spectrum. 



6i, 9/, (9^5/, d'^Qi, and As to measure c^^^^ to within standard errors of 10 percent. While the 
complexity of the first four fields go linearly in the number of particles, Ag is more expensive. 
The fact that we can achieve such consistency with such a small number of particles is not 
only remarkable, but numerically quite convenient. We describe the details of the analysis in 
App. D.l, and here simply provide a summary of the results. 

The spatial dependence of measured c^omb plotted for A = 1/3 h Mpc^^ and A = 
1/6 /i Mpc~^ in Fig. 5. These regions were chosen to maximize numerical stability as described 
in App. D. Fitting a constant to the value of c^omb displayed values after the UV cutofi^ 

gives in units of c^, we obtain 

cL^b(A = 1/3)= 0.96 ±.lxlO-^ (63) 
cL„b(A = 1/6) = 1.26 ±.lxlO-^ 

The RG flow between the two measured values is consistent with the prediction from 
perturbation theory as seen in Fig. 3. Furthermore the measured value from Consuelo agrees 
nicely with that predicted from matching to the nonlinear CAMB power spectrum. 
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Dependence of £'"i-i.iinb ^e\i for A— 1/3 (h/Mpc), A:,en— -16 




0.002 0.005 0.010 0.020 0.050 0.100 
^ext (h Mpc-i) 



Fi gure 4; In this plot we present to values obtained foT c^^^^^ as a function of the external momentum 
used in (62). We see that only as /cgxt 0, c'^ojnb becomes /Cext independent. This is so because at 
high k any higher derivative terms suppressed by powers of k/k]\jL are important. 
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Fi gure 5: Measurement of c^qj^jj in the UV with A = 1/3 /i Mpc^^ and A = 1/6 /i Mpc~^. 



5 Results 

The 0{6f) result of the computation of the power spectrum with our EFT is presented in 
Fig. 6. On top, we plot the ratio of the one loop power spectrum compared with the non-linear 
fit provided by the CAMB software with high precision settings, evaluated with the following 
cosmological parameters: Q\ = 0.75 ,Qm = 0.25, Qb = 0.04, h = 0.7, Ug = 1. The linear 
power spectrum is also obtained from CAMB with high precision settings. We take these 
data for all perturbative calculations done in the paper. Often in the literature the results of 
perturbation theory are plotted as ratio of the perturbation theory result versus a no-wiggle 
power spectrum. In this way the oscillatory features are still present in the plot, though they 
come mostly from the linear theory. We give this in the bottom part of the plot. 



24 



1.10 



1.05 



z 



0. 1.00 



^3 



0.95 



0.90 



Comparison with nonlinear power spectrum (CAMB): 
nonlinear power spectrum normalized. 



-i — Linear 

.* — Non-linear (CAMB) 

-t EFT (k;ren=- 16 h Mpc-^) 




« — • — • — • — • — ♦ 



0.05 0.10 0.15 0.20 
k (h/Mpc) 

Comparison with Nonlinear Power Spectrum (CAMB): 
"No-wiggle" normalized. 



0.25 0.30 




0.05 0.10 0.15 0.20 
k (h/Mpc) 



0.25 0.30 



Figure 6: The order 6f prediction from our EFT is compared with the CAMB non-linear output 
in the top, and to the no-wiggle power spectrum in the bottom, as well with the linear theory and 
Standard Perturbation Theory (SPT). The results from the EFT agree at percent level with the non- 
linear theory up to k 0.24/i Mpc~^ , when some high scale power seems to be missing. Results 
should improve already by going to 5^ order. The results are remarkably better than using SPT. The 
no-wiggle power spectrum we use is given by P5<5,no- Wiggle = 5.1 • 10^glog^(13g' + 2e)/(54 (7^(14 + 
731/(457g + 1)) + log(13g + 2e))2. 
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This plot is obtained after renormalizing the EFT prediction to match the power spectrum 
at kren. — 0.16/i Mpc~^. In detail, we perform the calculation at several increasing values of 
A, we choose the value of c^^^^^lao, A) to match the simulations' non-linear output, and then 
we extrapolate to A — )■ oo. This gives us c^^^^^^^^aQ, A = oo) ~ 6.2 x 10^^. Notice that naive 
considerations of virialization gave an estimated value for c^comh of order 10^^ [4]. The obtained 
smaller numerical value fits well with the decrease in the transfer functions for wavenumbers 
that are higher than the equality scale. 

The result for the power spectrum agrees at percent level with the CAMB non-linear fit 
up to A; ~ 0.24/i Mpc~^, where the EFT prediction begins to be smaller than the A^-body sim- 
ulation result. Results obtained with the approximate time dependence described in App. C 
are close to these ones, at percent level. As discussed in App. C, it is potentially dangerous to 
trust this approximation at high /c's for percent level precision, and luckily it is not very hard 
at all to perform the correct perturbation theory. We expect that the inclusion of the stochas- 
tic pressure and of higher order diagrams should improve the fit in the UV, possibly allowing 
us to fit the simulations to even higher /c's. Notice how the counterterm P13 -2 decreases 

' comb 

the power spectrum, compensating for the overshooting of SPT. It is difficult to interpret the 
percent disagreement that we have at moderate slow scales such as at ~ 0.12/i Mpc~^. At 
face value, it looks like that the computed power spectrum presents oscillations that are too 
large. These disagreements might be improved with the inclusion of higher order terms, or 
it might even be that at this level of precision, the results from the CAMB non-linear fit or 
from the A'^-body simulations might not be precise enough. These same improvements should 
reduce also the dependence on the renormalization scale: by changing the renormalization 
scale our results change by about 2%, at A; ~ 0.24/i Mpc~^. This dependence can be taken as 
a measure of the contribution from higher order terms. 

6 Conclusions 

Large scale structure surveys have the potential of becoming the next leading observational 
window on the physics of the early universe, potentially greatly improving what we are already 
learning from the CMB. Large Scale Structure physics is however much more complicated than 
the CMB due to the presence of large matter clustering at small scales. Since at non-hnear 
level different scales are coupled, these non-hnearities affect even large scale perturbations 
that are mildly non linear and so potentially treatable in a perturbative matter. In this 
paper we have developed the effective field theory of cosmological large scale structures in 
order to achieve a reliable predictability. Calculations in the effective theory are performed in 
k/k]\fL. The effective field theory is a cosmological fluid description for cold dark matter, and 
by extension all matter including baryons which trace the dark matter. The microphysical 
description is in terms of a classical gas of point particles, which we have smoothed at the 
level of the Boltzmann equation. We have exhibited and computed the various couplings that 
appear in the effective field theory, namely pressure and viscosity by matching to A"-body 
simulations, finding ~ 10~^ c^, etc. . We have developed the perturbative expansion for the 
power spectrum, which we have carried out at the 0{Sf). The fiuid parameters arise from 
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UV modes and alter standard perturbation theory. We have found that the corrections lead 
to a power spectrum in percent agreement with the full nonlinear spectrum as obtained by 
CAMB up to A; ~ 0.24h Mpc-\ 

It is a peculiar coincidence that k ~ 0.2Ah Mpc~^ is also the maximum k at which the 
popular technique of Renormalized Perturbation Theory (RPT) [9] works. While RPT is 
a very nice technique to compute non-linear corrections to the power spectrum, we stress 
that our approach is different at a qualitative and a quantitative level. At a qualitative 
level, RPT tries to solve, as exactly as possible, non-hnear equations for a pressureless ideal 
fluid. In our approach, instead, we try to solve non-linear equations for a different fluid. 
This has quantitative effects, as it is shown from the fact that as A — )■ oo our effective 
parameters like c^o^j^ do not vanish. What is more important about our EFT is that in can be 
improved. By performing higher order computations and by adding suitable counterterms, in 
principle arbitrary precision for reconstructing the power spectrum, or indeed any dark matter 
observables, can be achieved by going to a sufficiently high order in perturbation theory, 
on scales k < /catl- Techniques such as RPT or the renormahzation group approach [15] 
for example are still very nice techniques to perturbatively solve some non- linear equations, 
resumming many diagrams. It would be interesting to apply those techniques to solve the 
equations of our EFT. We leave this to future work. 

The effective field theory approach to large scale structure formation is complimentary 
to A^-body simulations by providing an elegant fiuid description. This provides intuition for 
various nonlinear effects, as well as providing computational efficiency, since the numerics 
required to measure the fluid parameters are expected to be computationally less expensive 
than a full scale simulation. Indeed we are able to achieve excellent agreement with the 
small down-sample of the full simulation we examined here. Of course, since the couplings 
are UV sensitive, it still requires the use of some form of A^-body simulation to fix the 
physical parameters, either by matching to the stress-tensor directly or to observables. But 
this matching is only for a small number of physical parameters at some scale and then the 
constructed field theory is predictive at other scales. 

There are several possible extensions of this work. A first extension is to go beyond the 
one-loop order to two-loop, or higher. This will require the measurement of several new 
parameters that will enter the effective stress-tensor at higher order, including its stochastic 
terms. Another extension is to compute the velocity fields and to include the small but 
finite contributions from vorticity, or to compute higher order A^-point functions, which can 
probe non-Gaussianity. Finally, another extension is to consider different cosmologies; in this 
work we have presented results on dark energy in the form of a cosmological constant. But 
one could equally consider other models for dark energy. This would presumably alter the 
value of the fluid parameters in a way that could be measured either from new simulations 
or determined by observation. Hopefully, our effective field theory for large scale striictures 
will help us use large scale structure surveys to uncover the physics of the beginning of the 
universe. 
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Appendix 

A UV Stress Tensor directly from Short Modes 

It is useful to write down the expression (22) directly in terms of short wavelength fluctuations. 
In this Appendix we provide these expressions. We define 

(64) 

(65) 

(66) 
(67) 

(68) 

where pj(x) = mai';(x). Note that cr*-' ^ cr'-' — a'^ ^ but they are related as follows 

= [^:^]a + {Pmv\v\\^ + h^(7r^- - prr^vl) + ^/(tt^ - p^v^{)\ ^ . (69) 

The second term is approximately piv\v\ (so it approximately cancels with —piv\v\ in kJ-') 
and the final term is small (as it is an overlap between short and long modes). Following the 
methods of [4] we obtain 



o1 


= m ^ J 


d^p(p^-pKx))(p^-p?(x))/(x,p) 






-^KXn))«-t'/(Xn))<5gHx-X„; 


(t>s,n 

di(f)s 


n 

= (t>n- (t>l,n 

= di(t)s,n 
n 






= di(t)s dj(f)s 





= [^y].+^^+o(ij) . (70) 
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Similarly, one can prove that satisfies 



+ 



dmdk(pldmdk(l)lS'^ - 2dmdi(j)ldmdj(pi 



So altogether we obtain the effective stress-tensor 



+ 



1 



where 



A2 



+ 



+ 



1 



(71) 

(72) 

(73) 
(74) 



We see that [t^^a is sourced by short wavelength fluctuations plus higher derivative correc- 
tions. 

Note that by taking the derivative dj this leading piece becomes 



with 



[psdi(f)s]A = ^di(f)s,n{^n)WA{x - X„) - [pidi(j)s]A , 

where the first term in (76) is given by 



(75) 
(76) 



^ m a ^9i05,^i(x„)iyA(x - x„) 



E 



vn?G [Xn — XfiY 



Erfc 



Alx, 



n -^n 



V2 



+ 



47r|x 



A2 



WAi^n - x^) 1^a(x - xO , (77) 



and the second term in (76) can be expanded as 

1 



[Pldi(f)s]A 



2A2 



(78) 



and this term should be included since it involves the background piece pj,, and so it includes 
a first order contribution. 
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B Expression for (5*^^^ 

The iterative solution for 5^^^ is given by 

- 256^ (^^) 

4Pi(a) (^J £q J d^p/3{p,q- p)6si{q- p)6si{p)P{q,k - q)6si{k - q) 

J A j d'pP{p,k - q- p)5s,{k - q- p)5sMmk - q)5s,{q)^ + 
2V2{a) J d^q j d^p/3{p,k-q-^5s,{k-q-p)6s,{^a{q,k-^6s,{q) + 
153(0) J '^^Q J d^pa{p,k — q — p)dsi{k — q — p)dsi{p)a{q,k — q)dsi{q) + 
1^4(0) J J d'^pa^Pjq — p)Ssi{q — p)Ssi{p)a{q,k — q)Ssi{k — q) + 
2V^{a) d^q J d^pa{p,q- ^Ssi{q- p)Ssi{^mk - q)Ss,{k - q) 

J d^q J d^pa{p,k - q- p)dsi{k - q- p)5si{p)P{q,k - q)5si{q)^ + 

2'De{a) J d^q J d^pP{p,q- p)Ssi{q- p)Ssi{p)a{q,k - q)Ssi{k - q) , 

where "Dj's represent the resuh of the integration of the Green's functions and the other 
time-dependent coefficients. They are given by 

pa pa 

Vi= ddd^G{a,h)n^{d)D\h) j ddd^V?{d)daG{d,d)D\df , (80) 

V2 = Snl^rn r dd r da a^'H\d)D'{dfD{d) G{d, d) + 2V^ , 

Jo ^ Jo 

i>3= r da 



dd I [3HlnmD{df + 2d^H'^{d)D'{df] [3nln^D{d)G{d,d) + 2d^H^{d)D'{d)daG{d,d)] , 



Vi = V^-A j ddd^ G{a, d)n\d)D{d)D'{d,f , 
Jo 

V5 = 3nlVlm r ddd^G{a,d)H\d)D'{d) f dd D{df^^^^p^ 
Jo Jo ^ 

-2 [ ddd^ G{a,d)H'^{d)D{d)D'{df + 2Vi , 
Jo 

Vg^V2 + 2 [ dda^ G{a, d)H'^{d)D{d)D'{df . 
Jo 



30 



Notice again the great simplification that occurs due to the fact that the growth factors and 
the Green's function do not depend on k, so that the time integrals and the momentum 
integrals decouple. 



C 0{Sf) Power Spectrum with Approximate Treatment 

In the main part of the paper, we performed perturbation theory with our EFT in a rigorous 
and exact way. However, it is possible to perform an approximate treatment that simplifies 
quite a bit the formulas. In this way it is simpler to follow the derivation and we therefore 
present it here. 

If the universe were to be EdS, then the solution for ^^^^ would be 5^"^ oc a" ~ D(a)". 
Thanks to this, all formulas simplify remarkably. Our universe is of course not of the EdS 

form, because of the cosmological constant. But it is tempting to extend the results obtained 
in EdS to the ones in ACDM universe by replacing in the EdS formulas the EdS growth factor 
with the growth factor in ACDM. In this we obtain 

{Si{k, ao)Siiq, ao))i_ioop - {27r)' S^'\k + q) (p22{k) + P^^{k) + ^^^(A;)) . (81) 
Here P22 and P13 are the time-independent one loop contributions given by 

/•OO / q r + 1 12 \ 

/ ^ (^^ - i^^^ + 2) log 42r^ + lOOr^ + ^ - 158 Pii,i{kr, A) . 

Jo 1 — r r J ' 

The counterterm contribution is given by 

A3..^..W-- f^[^^i:f '^-^ufeA). (83) 



where the time dependence of c^comb inferred in perturbation theory from (30) to be 



Notice that, within this approximation, the k dependence and the time dependence of c^comb^i^^ 

(3) 

is the same as the one of the source of 61 in the high k limit. 

This approximation is quite a good numerical approximation, and here below in Fig. 7 we 
present results of comparisons for P22 and P13, where we see that the disagreement is at percent 
level. This ratio is plotted in Fig. 8. Notice that at A; ~ 0.24 Mpc~^ where the non-linear 
corrections are of order of a few ten percents, a percent error in the loop calculation leads to 
an error that is dangerously close to one percent. Therefore, for percent precision, using this 
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approximate treatment at the high fc's that we can reach with the EFT corresponds to pushing 
the boundaries of safety. It is further important to stress that this is not a parametrically good 
approximation in 6i, as is our loop expansion. As noted in [29, 30, 31], this is an expansion in 
the smallness of the ratio ilrnia-)/{dlogD/d\ogay. A good fit is {d\og D/d\ogay ~ fi(a)^-^, 
explaining the possibility to make this approximation. Luckily, we find it not so hard to 
implement directly the correct perturbation theory, which is the one we present in this paper. 
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Figure 7: We plot the ratio as a function of time of P13 and P22 ols obtained using the rigorous 
time dependence and the approximate one. The ratio goes to 1 at early times when dark energy is 
irrelevant and the approximate treatment becomes exact. The results at redshift zero agree well, at 
percent level. 
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Figure 8: Plot of the ratio ^lm{o) / {dlogD / dloga)"^ . This being equal to one would justify the 
approximations done in this appendix. We see that at late times the difference is quite large. Since 
this is a correction to the one-loop term, this error can be acceptable for a one-loop calculation. On 
the contrary the approximation can become more harmful if one goes to higher loops. The actual 
result on the power spectrum of the approximation is even better than what shown in the plot, as 
more of the clustering happens before dark energy domination. 
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D Measuring in A/^-body simulations 



D.l Efficient calculation 



We break the calculation into three parts: the calculation of primary fields, the calculation 
of secondary fields, and the calculation of correlations. Primary fields are dependent solely 
upon the positions or velocities of the simulation particles, and secondary fields are dependent 
upon the primary fields, and we care about correlations between particular secondary fields. 
The immediate goal is to calculate the following expressions: 



clir) 



PAod^Pse ~ PAsd^Pee 
id^Pse)^ - d^Pssd'Pee 

^55 



PAsd^Pse - PasO^Ps. 
id^Pse)^ - d^Pssd^Pee 



(85) 
(86) 



based upon the two-point correlation functions: Pij{r). These two point correlation functions 
can be seen as the expectation value of the product of field / with field J, but with I evaluated 
at all points a distance r from all points r': 



Pijir) = {Iir + f')Jif'))air',r),r'= lim rrr^ r dr' dQ' dn/' Iif+ f') Jif'). 



(87) 



In the quantities being calculated the overall normalization cancels. The cost of evaluating 
each of these fields at all positions is prohibitive, so we approximate these correlations by 
measuring some large N number of pairs of points (pseudo)-randomly chosen but at fixed r: 

Pijir)^N-' J] Iifa')Jin') (88) 

{ra',rb'}eNset{r) 

where A^sct('") is a set of N pairs of points randomly selected from the space to be separated 
by a distance r. 

The following are the specific correlation functions used: 



PAsir 
PAoir 
Pssir 
Pseir 
Pe&ir 

'S5 

'seir 
ee(r 



d^P.sir 



{Asir + f')Siif'))a(r',r),r' 

{Mf+f')eiir'))^^r',r),r' 

{Siif+f')6iif'))Q^r',r),r' 

{Siif+f')eiif'))n(^r',r)y 

{eiir+r')eiif'))air',r)y 

{d^6lir + f') Slif'))n{r',r),r' 

{d^6iir + f') Qiir'))n{r',r),r' 

{d^eiif+f')eiif'))a(r',r),r' 



(89) 
(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 



These are functions of the gravitational short-mode field Ag, the over density S, and the 
velocity divergence 6, and relevant spatial derivatives. These fields are to be calculated from 
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the observed (or simulated) positions and velocities of point-sources at a fixed moment in 
time (redsliift=0 initially) . We smear the observation of the position of these sources with a 

gaussian, introducing a (soft) ultraviolet cutoff. 

To calculate cj, and in terms of Ag, d, 0, d'^d, d^Q we use a number of secondary fields. 
While it would be possible to numerically estimate the relative necessary spatial derivatives, 
it is simple enough to explicitly carry the operations out analytically, and treat them as 
independent secondary fields, avoiding the introduction of specious numerical error. 

(97) 
(98) 







[e](r) = 


1 ^ 

i=l 








i=l 


[9=^9] (f) = 


i=l,j=l 


[didj(l)s]{r) = 





(99) 
(100) 

(101) 
(102) 



[didjKij] (r) = [didjaij] (r) - [djin] (r) [diVj] (r) - [dim] (r) [djVj] (r) (103) 

- [7ri]{r) [didjVj]{r) - [didj7ri\{r) [vj]{r) 

[vm = N(r1/[p](r-) (104) 

[d,v,]ir) = [a,7r,](f)/[p](f) - h-](f)[a,p](f)/([p](f))^ (105) 
mvk]{r) = [d,d,7r,]{f)/[p]{f) ~ M(f)[9,p]/([p](f))2 - [d,7r,]{fmp]/{[p]{f)y (106) 

- M(f)[9.9,p](f)/([p](f))2 + 2[7r.](f)[9.p](f)[a,p](f)/([p](f))3 

mv,]{r) = ([p](f))-^{ - 6 {[dip]{r)r [d,p]{r) [7:j]{r) (107) 

+ 2 [mp]{r)r [djnj]{r) + [dApKr) [djp]{r) [nj]{f^ 

+ 2 [p](f) [9,p](f) ([a,7r,](f) [d,p]{f) + [d,d,p]{r) [7r,](f))) 

- {[p]{r)r (2 mp]{r) [9,7r,](r) + 2 [9,9,7r,](r) (M(r-)) 

+ [9.9.p](f) [9,7r,](f) + m7r,]{r) [d,p]{r) + [5^a,p](f) [7r,](f)^ 
+ ([p](rl)^[a,a,a,7r,](rl} 

These are, in turn, defined in terms of the following primary fields 
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(109) 


[didjp] (f) 
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n(iri{r) 
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[didjaij] (r) = 


mK'/{a\2T:fl^) J] didjW{r - rn){vnUvn)j 

n(iri(r) 


(116) 


[di[pmdj(f)]As] (r) = 




(117) 



where /i = r^ax and r]{'r) — {n s.t. |f — fn\ < Tmax}, and we have introduced the following 
functions for notational convenience: 



W{r) = exp J (118) 

diW(r) = -A\W(r) (119) 

didjW{r) = A^W{r) {A\rj - 5'^) (120) 

{di)%W(r) = A^ {diW{r) {A\rj - 5'^) + W{r)A\h rj + ^V,))) (121) 

i^n)j = X] (^n - rm)j r{fn - fm) [ ^ ^ ^ erfc(A|f; - fm\/V2) 

+ \ -AW{rn-rS\ , (122) 



where r)n — {m s.t. \rm — r„| < rmax}- We set rmax = 7A. It is efficient to precalculate and 
store {Cn)j for a given downsample of particles and A. Given that each simulation particle 
has a limited radius of influence rmax, a variety of parallelization strategies are available. We 
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Figure 9: The numerator and denominators of c^omb ^ measured with smoothing parameter 
A = 1/3 (left) and A = 1/6 (right), scaled to similar heights. This allows us to choose a con- 
venient region of measurement to avoid zero over zero contamination. Precision calculations 
in the future should extend measurements farther into the IR. 



have presented these fields, in detail, to emphasize that all such secondary fields rely on a 
relatively small number of primary fields, and it is only the primary fields which need concern 
themselves with the explicit number of particles down sampled from the simulation. 



D.2 Stability of Measurement and Region Selection 

The fiuid parameter of interest c^omb calculated as a ratio of polynomial functions of corre- 
lations. In fig. 9 we plot the numerator and denominator of these ratios for A = 1/3 [h/Mpc) 
and A = 1/6 {h/Mpc). It is worth noting strong confirmation of the effective field theory 
description is how well the numerator and denominator of the c^omb ratio tracks each other. 
With these sorts of statistical measurements, one should be careful of small fluctuations caus- 
ing misleading signal near zero over zero regions. For the measurements described in the 
paper we select a region where the denominator is 3a above zero. 



E Renormalizing at finite A 

The procedure as outlined in the main text involves renormalization for A = oo at some 
chosen k^en by fitting P}s~^°°^ for c^omb(^) against observation at that particular fcren- The 
power spectrum at all other k become predictions of the EFT. We could imagine to perform 
the same procedure at finite A. In this case, however, higher derivative terms suppressed by 
powers of k/A should be included. These terms do indeed vanish as A — )■ oo, but at finite A 
and k they are not negligible. In fact, as shown in Fig. 10, without the addition of higher 
derivative terms, the power spectrum deviates from the A = oo as — )■ A. One can indeed 
check that the values of c^^^-^ that is obtained by fitting in this way is off with respect to the 
correct value as obtained from running down to finite A the value of c^o^nb at A = oo by about 
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15%, depending on the cutoff used, see Fig. 11. Instead, by allowing for higher derivative 
terms, the correct value of c^omb is derived. 



A dependence of PsiSi (k,A) at^ren=0.16 




0.00 0.05 0.10 0.15 0.20 0.25 

k(h Mpc^) 

Figure 10: Prediction of the non-linear power spectrum without the addition of higher derivative 
terms, as we send A — )• oo, normalized to the non-linear power spectrum. We see that if we keep A 
finite, non-included higher derivative terms that scale as powers of k/A are important. Indeed the 
results improves as A = oo, which is the correct procedure. 
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Ratio of c^comb from direct fit at to c^comb from running down from A=ca at kj 

1.15 




^ren (h/Mpc) 

Figure 11: Ratio of the value of cl^^^{A,kren) as obtained from TUTminCI from C^omb('^ — ^ren) 
versus the one obtained by fitting directly the result of the EFT at finite A without the inclusion of 
higher derivative terms in k/A. At low k^en, the error is particularly pronounced for A = 1/6, as in 
that case /cren/A is not very small. Inclusion of higher derivative terms reduces the mismatch to few 
percent. 
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